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I. INTRODUCTION 

This paper continues the study of localised time- 
periodic solutions of the parametrically driven damped 
nonlinear Schrodinger equation, 

i^t + ^xx + 2|V^pV^ - = hilj* - i-ftp. (1) 

(Here /i, 7 > 0). Equation ([T]) is an archetypal equation 
for small and slowly- varying amplitudes of waves and pat- 
terns in spatially-distributed parametrically driven sys- 
tems. It was employed to model intrinsic localised modes 
in coupled microelectromechanical and nanoelectrome- 
chanical resonators [2 , solitons in dual-core nonlinear op- 
tical fibers [3 and dissipative structures in optical para- 
metric oscillators [4 . The discrete version of was 
studied as a prototype for the energy localisation in non- 
linear lattices [5 . (More contexts are listed in [1 .) 

In the previous publication [1 , its authors proposed to 
obtain the time-periodic solitons as solutions of the two- 
dimensional boundary- value problem with the boundary 
conditions 

ip{x^t) ^ as X ^ ±00; ip{x^t -\- T) = ip{x^t). (2) 

In the present paper, we apply this approach to the anal- 
ysis of complexes of solitons. 

Complexes (also known as molecules) are stationary 
or oscillatory associations of two or more solitons; they 
can be stable or unstable. Stable solitonic complexes, or 
bound states, were detected in a variety of soliton-bearing 
partial differential equations [6-15 . One mechanism of 
complex formation is the trapping of the soliton in a po- 
tential well formed by the undulating tail of its partner 
P", 8 . This mechanism is not accessible [7 to the para- 
metrically driven damped solitons though, as their tails 
are decaying monotonically . The exchange of resonant 
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radiation can also serve as a binding formula in nondissi- 
pative systems [Tj iQj , but in the damped-driven equation 
([1]) the radiation is nonresonant. A different mechanism 
was shown to operate here, which relies on the phase- 
stimulated growth or decay of the soliton's mass [10 . 

Bound states serve as long-term attractors in situations 
where there is more than one soliton present in the initial 
condition. For example, two like-polarity surface solitons 
in a vertically-driven water tank attract each other and 
form a stable bound state [11 . Unstable complexes do 
not have the same experimental visibility and can appear 
only as transients in numerical simulations. However, 
unstable complexes have a mathematical role to play: 
they work as the phase-space organisers [12 . 

The formation of complexes with an increasing num- 
ber of elementary constituents [13 gives rise to a higher 
degree of spatial complexity in the system, in the same 
way as the binding of shorter molecules into longer ones 
produces chemical compounds with increasingly complex 
properties. Previous analyses were confined to station- 
ary [10 and steadily moving [Mj [15] associations of the 
parametrically driven solitons. In the present paper we 
extend these studies to time-periodic complexes, thereby 
increasing the temporal complexity of the localised struc- 
tures. 

We consider time-periodic complexes as "stationary" 
solutions of Eq.Q on a two-dimensional domain — 00 < 
X < 00, < t < T. This allows to determine both stable 
and unstable complexes. Solutions of the boundary- value 
problem ([T]), ^ are path-followed in the parameter space 
— in the same way as free-standing periodic solitons were 
continued in the previous publication [1 . 

An outline of this paper is as follows. In the next 
section we describe bifurcations of the static two-soliton 
complexes. Of particular importance here are the Hopf 
bifurcations; these give birth to time-periodic solutions. 
We establish that the values of the damping coefficient 
are divided into two ranges. Namely, for 7 larger than a 
certain threshold, the complex suffers one or more Hopf 
bifurcations as h is varied. Below the threshold 7, no 
Hopf bifurcations occur. 



In section III, the Hopf-bifurcation points of the sta- 
tionary complexes are exploited as the starting points of 
the T{h) curves for the time-periodic complexes. These 
curves are traced as we continue the periodic bound 
states in a parameter. Depending on the number of the 
Hopf bifurcations suffered by the static complex, we have 
one, two or more branches of the periodic solutions ema- 
nating out of it. Complexes resulting from different Hopf 
bifurcations follow different transformation routes. 



In the concluding section (section IV) the results on 
stationary and time-periodic complexes are summarised 
in the form of a two-soliton attr actor chart. Included in 
this chart are also some quasiperiodic attractors. 



II. STATIONARY TWO-SOLITON COMPLEXES 

The two free-standing stationary soliton solutions to 
Eq.Q are distinguished by the subscripts + and — : 

iIj±{x) = A± ex.p{—i6±) sech(A±x), 



where 



,1.7 

9+ = -arcsm-, 



and 6- = 7r/2 — 6>+. The ip- soliton is unstable for all h 
and 7. The soliton ?/^+ is stable when the difference h — j 
is small but loses its stability to a time-periodic soliton 
when h exceeds a certain limit /iHopf(7)- 

The solitons ?/^+ and ip- can form a variety of bound 
states, or complexes [TOl [TH [15] . (For example, in the 
previous paper [1 we mentioned a complex ?/^(_+_), that 
is, a symmetric stationary association of two solitons ip- 
and one All complexes involving the ?/^_ solutions 

are, expectably, unstable; however two ip-^ solitons can 
form a stable bound state [10 . 

Previously, the two-soliton complex was known 

to exist only for sufficiently large values of damping [lOl 
[T5] . We have now established that this complex exists 
for all 7 > 1.5 X 10~^. Its domain of existence on the 
(7, /z)-plane is not bounded from above except that for h 
greater than 

^cont \/l +7^ 

the complex is unstable to the continuous-spectrum per- 
turbations (as any other solution decaying to zero at the 
infinities). Reducing h from /icont for the fixed 7, we ob- 
tain one of two possible types of bifurcation diagrams on 
the (/i, E) plane, where the energy E is defined by 



E ■ 



*2-| 



dx. (4) 



(The energy is not an integral of motion when 7 7^ 0; 
however, E is obviously a constant for time-independent 
solutions and can be used as a physically meaningful bi- 
furcation measure.) 
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FIG. 1. Energy of the stationary two-soliton complex and 
stationary multisoliton solutions obtained from this complex 
by continuation in h for the fixed 7. (a) 7 = 0.01; (b) 7 = 
0.4. Solid curves show stable and the dashed ones unstable 
solutions. Note two intervals of stability of the complex 
in (b). 



The diagram of the first type [FigjTJa)] arises when h 
is decreased for a fixed small 7 (7 < 0.292). In this 
case, there is only one turning point, h = /ign, with 
^sn = ^sn(7) > 7- [For the parameter value 7 = 0.01 
which we used to create Figjlja), hsn = 0.02972; for 
7 = 0.1, the turning point is at /ign = 0.25, and for 
7 = 0.25, hsn = 0.49.] As h approaches h^n along the top 
branch, the two-soliton solution develops a third 

hump halfway between the two humps that are already 
there, with the distance between the lateral humps re- 
maining unaffected by this development. The complex 
obtained by the continuation of this solution to the bot- 
tom branch can be identified as a three-soliton bound 
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FIG. 2. (Color online) The existence and stability domains 
of the stationary two-soliton complex The region of 

existence of the complex is bounded from below by the dotted 
line; in the upwards direction it extends beyond the value 
/icont, without bound. The complex is stable only in a small 
part of this region [tinted yellow (light gray)]. The stability 
domain is bounded by the curve h — /icont (7) on the top and 
by the lines of four Hopf bifurcations (1, 2, 3, and 4) on other 
sides. The dashed curve is the line of the Hopf bifurcation 
of the single soliton '0+. (The soliton is unstable above this 
line.) 



state ?/^( \ ). As we continue away from /ign along the 

bottom branch in Figjlja), the solitons bound in this 
complex (the two side solitons) diverge to the infinities 
on the X axis. 

All branches in the diagram of the first type consist of 
unstable solutions. (Our approach to the stability anal- 
ysis of stationary solutions has been outlined in [1 .) 

A somewhat different diagram arises for larger values 
of 7 (7 > 0.292), see Figjljb). This bifurcation diagram 
has been described in [10 for a particular 7 (7 = 0.565); 
here, we reproduce it for a different value of the damping 
coefficient. Reducing h from /icont for the fixed 7, the 
branch resulting from the two-soliton solution de- 
velops two turning points instead of one. As we pass the 
first turning point, the complex transforms into 

the solution. Moving away from the point along 

the bottom branch, the ?/^( ) complex acquires a third 

hump. This branch does not continue all the way to /icont 
but turns again, into a branch with even a lower energy. 
On this branch, the three-hump solution can be identi- 
fied as ^_). The lowest branch continues to the point 

h = J. This point defines the lower boundary of the 
domain of existence of the stationary complexes which 
result from the path-following of the two-soliton solution 
?/^(++). As we approach the point h = the distance be- 
tween the two side solitons in the ?/^( \ ) complex tends 

to infinity. 

For 7 just above 0.292, all solution branches in the 



diagram of the second type are unstable. However, as 
7 exceeds 0.34, a stability window opens in the '^(++) 
branch. The existence and stability domains of the two- 
soliton complexes on the (7, h) plane are shown in Figj2j 

We were not able to obtain a symmetric two- or three- 
soliton complex for 7 = 0. If we fix h and continue in 
7 towards 7 = 0, the separation distance between the 
solitons in the complex grows without bounds; hence we 
conjecture that symmetric multisoliton complexes do not 
exist for 7 = 0. (There are nonsymmetric complexes with 
7 = though; see [H].) 

The shape of the E{h) curve corresponding to 7 = 0.4 
[Figjljb)] looks similar to that of the E{h) curve for 
7 = 0.565 [1^. The main difference between the dia- 
grams pertaining to these two values of 7 is that when 
7 = 0.565, the stability region of the two-soliton solu- 
tion is seamless, i.e. does not have instability gaps in it, 
whereas in the 7 = 0.4 case, the stability region consists 
of two segments of the curve separated by an interval of 
instability. This difference is reflected by the shape of 
the stability domain on the (7, /i)-plane (note the "insta- 
bility bay" on the north-west coast of the stable region 
in Figj2|. 

Each point of the "coastline" of the stability "conti- 
nent" in Fig|2] corresponds to a Hopf bifurcation of the 
stationary complex [except for points along the curve 
h — /i-contlT)]- The "coastline" consists of four segments 
(marked 1, 2, 3, and 4 in Figj2|. Continuing in h along 
a vertical line 7 = const one crosses one, two or four of 
these; accordingly, for a given 7, the complex may un- 
dergo one, two or four Hopf bifurcations. 



III. TIME-PERIODIC COMPLEXES 

The first segment (marked 1 in Figj2]) is defined as 
the "south coast" of the tinted continent. It extends 
from 7 = 0.34 to larger 7 without a visible bound — 
presumably all the way to 7 = 00. The line of the second 
Hopf bifurcation (marked 2) represents the "north coast 
of the southern peninsula" in Figj2] it is bounded by 
7 = 0.34 on the left and 7 = 0.413 on the right. The 
"south coast of the northern peninsula" corresponds to 
the third Hopf bifurcation (marked 3); this extends from 
7 = 0.39 to 7 = 0.413. Finally, the top, fourth, Hopf 
bifurcation arises for 7 between 0.39 and 0.445 (marked 
4). When 7 is greater than 0.445, the complex undergoes 
just one Hopf bifurcation as h varies (the one marked 1). 



A. The first Hopf bifurcation 

In this subsection, we path-follow time-periodic com- 
plexes born in the lowest Hopf bifurcation (i.e. detaching 
from the south coast of the tinted "continent" in Figj2|. 
We take 7 = 0.565 as a representative value of damping 
in the region where the stationary two-soliton complex 
undergoes only one Hopf bifurcation and 7 = 0.35 and 
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0.38 in the region where there is more than one Hopf 
point. 

When 7 = 0.565, the (only) Hopf bifurcation is at 
hni = 0.94. Using this value as a starting point in our 
continuation process, results in the bifurcation diagram 
shown in Figjsj^a). In order to articulate details of the 
diagram, we supplement the graph of the period T{h) 
with a plot of the averaged energy, defined by 

E = ^I^E{t)dt, (5) 

where E{t) is given by Eq.Q. We also evaluate Floquet 
multipliers as described in jlji. 

At the starting point hni = 0.94, the Floquet spec- 
trum includes three unit eigenvalues and two complex- 
conjugate pairs with moduli smaller than one. As h is 
decreased from h = 0.94, two unit eigenvalues remain in 
the spectrum while the third one moves inside the unit 
circle along the real axis. This positive eigenvalue de- 
creases in modulus until it passes to the negative semi- 
axis at h = 0.92; once the eigenvalue has become nega- 
tive, it starts growing in absolute value. Eventually, as 
h reaches the value of 0.897, the negative real eigenvalue 
crosses through /i = — 1. A period-doubling bifurcation 
occurs at this point; as h drops below 0.897, the periodic 
complex becomes unstable but a stable double-periodic 
solution is born. Note that the destabilization occurs 
not at the turning point of the E{h) curve (which is at 
h = 0.89665) but for a slightly larger /i, i.e. before the 
turning point is reached. Figjsjb) shows a representative 
solution on the lower, "horizontal", branch of the E{h) 
curve. 

As for the two complex pairs, the eigenvalues consti- 
tuting one of these grow in absolute value as we move 
along the "horisontal" branch towards smaller h. At the 
same time, the imaginary parts of these eigenvalues de- 
crease and the pair converges on the positive real axis — 
just before crossing through the unit circle. The two real 
eigenvalues cross through /i = 1 almost simultaneously, 
as the curve turns back at = 0.89665; after that, they 
remain outside the unit circle. The other complex pair 
also converges on the real axis but remains inside the unit 
circle along the entire curve. 

As h is decreased and we approach the turning point 
in Figjsja), the amplitude of oscillations grows and the 
solution transforms into a sequence of soliton fusions and 
fissions. Two solitons merge into one entity which then 
breaks into two constituents, and this process continues 
periodically; see Figjsj^c). 

The whole of the "vertical" branch of the E{h) curve 
is unstable. The branch ends at the stationary ?/^(+++) 
solution (here = 0.901). As we approach the endpoint 
of this branch, the two real (positive) eigenvalues with 
/i < 1 and one of the two eigenvalues with ja > 1 move 
closer to 1. At the endpoint, the spectrum includes three 
unit eigenvalues and two real eigenvalues close to 1. This 
corresponds to the spectrum of a stationary three-soliton 
complex near its Hopf boifurcation. 



Proceeding to the region with more than one Hopf 
point, we consider 7 = 0.35. Here, the "lower" Hopf 
bifurcation occurs at hni = 0.805. This bifurcation is 
supercritical; as h drops below hni^ the stationary two- 
soliton bound state loses its stability to a periodic two- 
soliton complex which is born at this point. At the bifur- 
cation point, the spectrum of the Floquet multipliers in- 
cludes three unit eigenvalues and two complex-conjugate 
pairs inside the unit circle — one with Re /i < and 
the other one with Re /i > 0. As we continue the pe- 
riodic complex towards smaller /i, the negative-real-part 
pair converges on the real axis inside the unit circle, af- 
ter which one of the resulting negative eigenvalues grows 
in absolute value and, at = 0.79, crosses through 
/i = — 1. The periodic complex loses its stability to a 
double-periodic bound state of two solitons. As we con- 
tinue the unstable branch, it makes a number of turns 
(Figjij^a)), the spatiotemporal complexity of the solution 
increases (Fig|4]^b)) but it never regains its stability. 

Another representative value of 7 with two Hopf bi- 
furcations, is 0.38. Here, the continuation of the two- 
soliton complex from the lower Hopf point results in the 
T{h) curve similar to the 7 = 0.35 case (Fig|4]^a)). Like 
in the 7 = 0.35 case, the solution loses stability in a 
period-doubling bifurcation. We did not path-follow the 
unstable branch far beyond the bifurcation point. 

Summarising results of continuation from the first, 
"lowest" , Hopf bifurcation in Fig j2j we note that the bi- 
furcation is supercritical — both for large and small 7. 
Another common feature is the loss of stability resulting 
from a Floquet multiplier crossing through ja = —1. Since 
this bifurcation occurs before the first turn of the T{h) 
curve, it always gives rise to a stable double-periodic solu- 
tion. It is also fitting to note that all time-periodic com- 
plexes emerging in the first Hopf bifurcation are symmet- 
ric in space (i.e. invariant under the reflection x —x). 



B. The second and the third Hopf bifurcations 

When 7 lies between 0.34 and 0.39, the stationary two- 
soliton complex suffers two Hopf bifurcations, at hui 
and /iH2, with /ih2 > ^hi- (These are marked 1 and 
2 in Figj2j) In this subsection, we describe the continu- 
ation of periodic solutions detaching at /ih2 (the second 
of the two bifurcations) for several representative values 
of damping. 

The second Hopf bifurcation is subcritical: the emerg- 
ing periodic branch is unstable and coexists with the sta- 
ble stationary branch. That is, the periodic branch ini- 
tially extends down in /i, see the 7 = 0.35 and 7 = 0.38 
curves in Figjs] At some point the branch turns back 
after which h grows without any further U-turns; no- 
tably, it grows beyond the interval of the stable station- 
ary bound states (Figj5|. 

The periodic branch ends at an unstable stationary 
complex. The endpoint corresponds to the "concealed" 
Hopf bifurcation of the stationary solution where a pair 
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of complex- conjugate eigenvalues crosses from one half of 
the complex plane to the other but the solution remains 
unstable due to additional unstable eigenvalues. 

When 7 is set on 0.35, the whole periodic branch is 
unstable but for 7 as close as 0.36 a narrow stability 
window appears inside it. As 7 grows from 0.36, the sta- 
bility window expands — see the 7 = 0.38 curve in Fig|5] 
which features a sizeable stability interval hi < h < 
with hi = 0.9415 and /12 = 1.015. Within this stability 
window, the periodic complex has two complex-conjugate 
pairs of Floquet multipliers /i4 = /ig and /ie = with 
< 1 (in addition to two unit eigenvalues). As h is de- 
creased below hi, the first pair (/i3,4) moves outside the 
unit circle, with the second pair remaining inside; when 
h is raised above /i2, the unit circle is crossed by the sec- 
ond pair (//5,6)7 with the first pair remaining inside. Thus 
the stability interval is bounded by the Neimark-S acker 
bifurcation on each side. This observation suggests that 
a quasi-periodic two-soliton complex should be born on 
the crossing of either stability boundary, hi and /12 — 
the conclusion confirmed by direct numerical simulations 
of Eq.Q. (Quasiperiodic solutions can obviously not be 
captured by the periodic boundary- value problem; the di- 
rect numerical simulation remains the only feasible way 
of determining them.) 

It is worth mentioning here that the periodic two- 
soliton complexes coexist with periodic one-soliton so- 
lutions. (For example, for 7 = 0.35, the periodic free- 
standing soliton exists between h = 0.75 and h = 1.02; 
see Fig. 2 in [1 .) However the one- and two-soliton 
branches are not connected. 

When 7 is between 0.39 and 0.413, the stationary com- 
plex undergoes four Hopf bifurcations, hm < /ih2 < 
hm < ^H4 (marked 1, 2, 3, 4 in Figj2j). This is the 
interval of 7 that contains the top "peninsula" in Fig|2j 
Choosing 7 = 0.41 as a representative value of damping, 
we path-followed the periodic complex which is bifurcat- 
ing off at the point /ih2- Like in the case 0.34 < 7 < 0.39 
discussed above, the bifurcation is subcritical: the emerg- 
ing periodic branch is unstable and initially extends down 
in h. As in the previously discussed case, the T{h) curve 
turns {ai h = htn = 0.9447) and the entire subsequent 
continuation proceeds in the direction of increasing h 
(Fig©. 

To describe the motion of the Floquet multipliers, it is 
convenient to start somewhere within the "upper" part 
of the 7 = 0.41 branch, e.g. at = 0.99. At this 
point, the spectrum of linearisation includes two pairs 
of complex-conjugate multipliers /i4 = /ig, /ie = both 
with < 1. However, in contrast to the previously dis- 
cussed scenario, neither of these two pairs crosses through 
the unit circle as h is inceased or decreased and so the 
periodic complex with this 7 does not experience any 
Neimark-S acker bifurcations. 

As h is decreased from 0.99, the multipliers /is ^4 con- 
verge on the real axis and, at the turning point /itn, cross 
through /i = 1 (almost simultaneously). The other com- 
plex pair, /is, 6 5 remains inside the unit circle. Therefore, 



the turning point corresponds to a saddle-node bifurca- 
tion of limit cycles. If we, instead, increase h starting at 
h = 0.99, it is the jh^^q pair that converges on the real 
axis, just before /is becoming equal to one. At this point 
the periodic branch rejoins the branch of stationary com- 
plexes; this value of h is nothing but /ih3, the point of 
the third Hopf bifurcation of the stationary two-soliton 
bound state. At h = hus^ the Floquet spectrum includes 
three unit eigenvalues, a real eigenvalue /ie close to (but 
smaller than) one, and a complex-conjugate pair /i3,4 in- 
side the unit circle. Thus the periodic complex remains 
stable in the whole range between the turning point htn 
and the point hus where it rejoins the (stable) stationary 
branch. 

Since the saddle-node bifurcation point htn lies below 
/iH2, there is an interval htn ^ h < /ih2 where we have 
bistability between the stationary and time-periodic two- 
soliton complexes. 

In summary, the second Hopf bifurcation is always sub- 
critical; the continuation connects it either to the third 
(supercritical) Hopf bifurcation, or to a concealed bi- 
furcation of unstable two-soliton complexes. All time- 
periodic solutions arising in these bifurcations are sym- 
metric in space. 

C. The fourth, symmetry-breaking, Hopf 
bifurcation 

The locus of the fourth Hopf bifurcation is a stretch of 
the north-west coast of the "continent" of stable station- 
ary complexes in Fig j2] (marked 4). The "north-western 
coastline" extends from 7 = 0.39 to the point 7 = 0.445, 
where it meets the continuous-spectrum instability curve 
h = /?^cont(7)- At the bifurcation point a pair of complex 
eigenvalues A, A* of the eigenvalue problem 

{n-jj)p{x) = xjp{x), 

cross through the imaginary axis. Here H — 7 J is the op- 
erator of linearisation about the stationary solution (see 
for details). The bifurcation is symmetry breaking: 
unlike three other Hopf bifurcations, the corresponding 
eigenfunctions p(x) and p*(x) are odd (antisymmetric): 
p{—x) = —p{x). Accordingly, the time-periodic solu- 
tions which are born in this bifurcation describe out-of- 
phase oscillations of two identical solitons making up the 
complex [see Figj6]^b)]. 

The fourth Hopf bifurcation is supercritical: the 
emerging periodic branch is stable and extends up in 
h. For 7 = 0.41 which we take as a representative 
value of damping, this bifurcation occurs at /ih4 = 1.037. 
At the bifurcation point, the Floquet spectrum includes 
three eigenvalues /ii,2,3 = 1, two real positive multipli- 
ers /i4^5 < 1 and several complex- conjugate pairs with 
|/i| < 1. As h grows from hn^, one of the unit eigenvalues 
moves inside the unit circle, but as h is further increased, 
it reverses and moves out. At this point {h = 1.049), 
a saddle-node bifurcation of cycles occurs; the periodic 
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solution loses its stability and the branch turns back 
[Figj6|a)]. As we continue it further, the real and com- 
plex eigenvalues move back and forth through the unit 
circle; some pairs converge on the real axis — but the 
solution never regains its stability. 

After a lengthy excursion into the (h^T) plane 
[FiglGja)], the periodic branch rejoins the branch of (un- 
stable) stationary complexes (at hcnc = 1.082). 
The spectrum of the endpoint stationary solution in- 
cludes three unit eigenvalues, a positive eigenvalue fi^ < 
1, and two complex-conjugate pairs, with |/i5,6| < 1 and 
|/i7^8| > 1- The value h = hcnc pertains to the con- 
cealed Hopf bifurcation of the unstable stationary com- 
plex 

IV. THE TWO-SOLITON ATTRACTOR CHART 
AND OPEN PROBLEMS 

Figjr] summarises our results on the stationary and 
periodic two-soliton attractors. This diagram comple- 
ments the single-soliton attractor chart compiled in the 
first part of this project |1 . The two-soliton chart is 
in qualitative agreement with results of direct numerical 
simulations [16 . 

In FigjTj we have included stable quasiperiodic 
complexes [highlighted in purple (dark gray)]. The 
boundaries between the stable-periodic and stable- 
quasiperiodic domains are defined by the Neimark-Sarker 
bifurcations of the periodic complexes; these admit an ac- 
curate demarcation using our method (i.e. by monitoring 
the Floquet multipliers along the periodic branches) . On 
the other hand, in order to determine where the stable 
quasiperiodic solution ceases to exist, we had to relin- 



quish our continuation approach in favour of direct nu- 
merical simulations of Eq.(p!]). We have performed only 
a few runs and hence Figj7[gives only a schematic posi- 
tion of the outer boundary of the quasiperiodic stability 
domain. In order to demarcate this boundary more accu- 
rately, one would need to perform numerical simulations 
more extensively. This is beyond the scope of our present 
study. 

The region of bistability of stationary and periodic 
complexes also needs to be accurately delimited. So far, 
we have only demarcated a small portion of it; see the 
black mark in FigjTj 

Finally, it would also be interesting to continue peri- 
odic solutions bifurcating from the stationary complexes 
in the "concealed" Hopf bifurcations, where the station- 
ary solution remains unstable on both sides of the bi- 
furcation due to additional eigenvalues with positive real 
parts. In our continuation process bifurcations of this 
sort would typically arise as the endpoints of the peri- 
odic branches starting at the proper Hopf bifurcations of 
the stationary complexes. Starting at the concealed bi- 
furcations would produce an additional wealth of periodic 
branches some of which may have stable segments. 
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FIG. 3. (Color online) (a) The average energy (main panel) 
and the period (inset) of the periodic solution arising for 7 = 
0.565. The solid curve shows the stable and the dashed one 
unstable branch. The point of the period-doubling bifurcation 
is marked by a vertical arrow, (b-c) Representative solutions 
on the lower branch in panel (a), (b): /i = 0.92, T = 13.973; 
shown is the absolute value of V^. (c): h — 0.90, T — 32.729; 
shown are the level curves of l-^l. The time interval covered 
by (b) and (c) includes three periods of oscillation. 




FIG. 4. (Color online) The first branch of the two-soliton 
time-periodic complex for 7 = 0.35 and 7 = 0.38. (a): the pe- 
riod of the solution. The solid curves show the stable and the 
dashed ones unstable branches. The circles mark the starting 
points of the continuation (the stationary complex 
(b): A two-soliton periodic solution with complex temporal 
behaviour arising at the end point of the 7 = 0.35 curve in 
(a). (Here h = 0.741, T = 15.9.) The figure shows just one 
period of oscillation. 
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FIG. 5. The second branch of the periodic two-sohton so- 
lution. The empty circles mark the starting points of the 
continuation — the points h = hu2 where the stationary two- 
soliton complex suffers the second Hopf bifurcation. The full 
circles mark the endpoints. The endpoint of the 7 = 0.41 
curve corresponds to the third Hopf bifurcation (/ihs); the 
endpoints of the 7 = 0.35 and 7 = 0.38 curves lie inside the 
instability domain of the stationary complex. The solid curves 
show the stable and the dashed ones unstable branches. 
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FIG. 6. (Color online) (a) The branch of periodic 2-soliton 
complexes oscillating out of phase with each other. The empty 
circle marks the starting point of the continuation — the point 
h = /iH4 where the stationary two-soliton complex suffers the 
symmetry-breaking Hopf bifurcation. The full circle marks 
the endpoint. The short solid segment near the beginning of 
the curve represents the stable solution; the rest of the branch 
(dashed) is unstable, (b) A representative solution on the 
stable part of this branch. Here h — 1.0493 and T — 1.991; 
several periods of oscillation are shown. 
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FIG. 7. (Color online) Two-soliton attractor chart. The chart 
is still under construction; the outer boundaries of the pur- 
ple (dark gray) area (populated by quasiperiodic two-soliton 
attractors) are still to be refined. The black mark (pointed 
to by the arrow) initiates the periodic/stationary bistability 
region. The dashed curve is the line of the Hopf bifurcation 
of the '0+ soliton, shown just for the reference purposes. 



